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We study the large mass asymptotics of the Dirac operator with a nondegenerate mass matrix 
m = diag(mi, rri2, ms) in the presence of scalar and pseudoscalar background fields taking values 
in the Lie algebra of the (7(3) group. The corresponding one-loop effective action is regularized by 
the Schwinger's proper-time technique. Using a well-known operator identity, we obtain a series 
representation for the heat kernel which differs from the standard proper-time expansion, if mi 
m,2 m3. After integrating over the proper-time we use a new algorithm to resum the series. The 
invariant coefficients which define the asymptotics of the effective action are calculated up to the 
fourth order and compared with the related Seeley-DeWitt coefficients for the particular case of a 
degenerate mass matrix with mi = 7712 = rriz- 
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The effective action plays the central role in lagrangian quantum field theory. Dividing fields into a classical 
background and quantum fluctuations and integrating out the last ones one obtains the effective action which properly 
accumulates the short distance dynamics of quantum fields In this letter we study the real part of the one- 

loop effective action with virtual heavy fermions of different masses using the Schwinger's proper-time technique 0. 
In QCD the nondegenerate mass matrix of heavy (constituent) quarks results from the spontaneous breakdown of 
chiral symmetry, as it is the case, for instance, in the Nambu - Jona-Lasinio model 0], or from the specially added 
' invariant term to the QCD Lagrangian that regulates the infrared behaviour of the effective action ||. We do not 
consider here the manifest chiral symmetry breaking effect on the effective action ||. A careful analysis of this 
' problem would lead us too far away from the subject, leaving the present result without changes. It might be well 
to emphasize, however, that a nondegenerate mass matrix of heavy fermions in theories with spontaneous breakdown 
of chiral symmetry appears as a consequence of manifest symmetry breaking by the nondegenerate mass matrix of 
light fermions. The mathematical formalism being presented in this letter is a necessary element of the approach 
which faithfully mirrors the vacuum structure of such a theory. This approach was formulated for the theory with an 
explicit and spontaneous breakdown of the global SU (2) x SU (2) chiral symmetry in || and generalized to all orders 
of the asymptotic expansion in Q). Here we extend this scheme to the SU(3) x SU(3) chiral theory where the mass 
matrix has a form m = diag(mi, ni2, m^). 

The explicit forms of dominant effective local vertices induced by virtual heavy fermions in general spontaneously 
broken gauge theories have been obtained in However both the method of our evaluation of the heat kernel as 

well as the result, which is cast as an asymptotic series for the effective action with the order by order chiral invariant 
structure of the derived asymptotic coefficients, are different from the cited papers and new. Let us clarify: we do 
. ' not use the equation of the Schrodinger type for the heat kernel in the first stage of calculations, prefering the direct 
expansion of the heat kernel in powers of background fields, like, for instance, it has been reviewed in fl(i|j . Additionally, 
' we suggest a new algorithm for resummations. Starting from this place our calculations essentially deviate from the 
known ones. The resummation procedure in j^] finally leads to an asymptotic series where the n-th term gives the 
complete 0{m~ n ) contributions (including all powers of lnm) to the effective action. In chiral theories this structure 
immediately breaks chiral symmetry at each order of n. The most direct way to see this is to consider the theory 
with a linear realization of chiral symmetry. In the spontaneously broken phase the symmetry transformations of 
scalar and pseudoscalar fields include mass dependent parts. Therefore, the asymptotic coefficients must have terms 
with different powers of masses, to be chiral invariant. The resummation procedure must take this fact into account 
properly. In our approach resummations are organized in such a way that every coefficient of the asymptotic series is 
automatically chiral invariant, if the Lagrangian possesses this symmetry. We consider this property to be a crucial 
condition on any generalization of the Schwinger - DeWitt result. Moreover, the resummation procedure used in Q is 
slightly misleading. Indeed, the terms like, for example, [<fi, m 2 ] where <fc is a background field, are considered to be of 
m 2 order. However, it is clear that the commutator contributes as the difference of mass squares, i.e., as ~ {m\ — m 2 ,) 
(in the case of 2 x 2 mass matrix m). This value can be small and as a result contribute to lower orders. 
A further approach consists in studying the proper-time asymptotics of the heat kernel with arbitrary matrix- 



*On leave from the Laboratory of Nuclear Problems, JINR, 141980 Dubna, Russia 



1 



valued scalar potentials and nontrivial mass matrices as a special case (see, for instance, |llL|l2[) with Seeley-DeWitt 
coefficients known up to n — 5 for an arbitrary operator of Laplace type. The general formalism and some references 
can be found in 13 1, chiral gauge theories are reviewed in |1C(| . One should not mix the subject of those investigations 
with the problem considered here. We do not study the proper-time asymptotics although in the limit mi = m-i = 7713 
our asymptotic coefficients coincide with the standard Seeley-DeWitt coefficients. In fact, we use the proper-time 
representation only to separate the field dependent part of the heat kernel from the mass dependent one. We do not 
expand the mass dependent piece of the heat kernel in powers of the proper-time; instead we manipulate all terms of 
the proper-time series at once. This is an essential part of our approach and it means that there is no way to obtain 
our result from the well known proper-time asymptotics. 

After these general remarks, to be more specific, let us consider an euclidean quantum field theory with a fermion 
propagator D^ 1 depending on the background fields collected in 7. In the one-loop approximation the real part of 
the corresponding effective action is given by 
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We use here the Schwinger's proper-time representation for the modulus of the functional fermion determinant ||. 
D^D is a second order elliptic differential operator, D^D = m 2 + B, with B — — d 2 + Y. The symmetry group acting on 
background fields is the SU (3) x SU(3) group of global chiral transformations. The mass matrix m = diag(mi , m 2 , m 3 ) 
can be written in the following way 



y^^mjEj, (Ei)jk = Sijdik, EiEj— SijEj. 



(2) 



We find it convenient to use the orthogonal basis, Ei, in our analysis. The integral over t is divergent at the low 
limit and needs to be regularized. This can be done by inserting the regularizing kernel p(t, A 2 ) with the ultraviolat 
cutoff parameter A. In the following we do not need an explicit form for p(t, A 2 ). Using the technique developed by 
Fujikawa 111] one can obtain 



W[Y] 



dt 
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where A — B — 2ipd/ \/t. As soon as we have derivatives it is necessary to clarify the meaning of the trace in Eq.(||). 
The coordinate space {x} is assumed to have no boundary, so that one can integrate by parts rendering the functional 
trace cyclic. Nevertheless, it is useful to define the cyclically symmetrized trace of matrix-valued functions Aj by 



stT(A 1 ,A 2 ,...,A n ) = V -tr (A1A2 . . . A n ) , 



(4) 



where perm means n possible cyclic permutations inside trace. For algebraic objects str is equivalent to tr provided 
the trace is cyclic. 

Since we do not want to expand the mass-dependent part of the heat kernel in powers of the proper-time, we need 
the operator identity, which is well-known in quantum mechanics 



tr e 



Here f n (t, A) is equal to 
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ds n A(s 1 )A(s 2 )...A(s n ), 
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where A(s) = e sm Ae s,n . For our purpose we need the first five terms of the series. The n-th coefficient is equal to 
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where Aj = EiA, and the totally symmetric coefficients c^^..,^ (t) are equal to 



c ijk 



t 2 \ A H A 
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with the definition Ay = to 2 — m 2 . It is not difficult to write down the expression for any coefficient Ci 1 i 2 ,..i n in (0). 
If the indices are all equal, one can obtain that Ci = c„ = c«...j. In the limit m\ = 7712 = 7713 = M this becomes 
c, = exp(— tM ) and the exponent is totally factorized. Thus one can obtain the standard inverse mass expansion. 
In the nondegenerate case for the first five terms in (0) we find 
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Substituting this expression in Eq.(|3|), putting there A = Y — d 2 — 2ipd/\/t and integrating over the four- momenta 
p M we obtain 
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Here we have used Y\ = EiY and considered only terms up to and including t 3 order (we count powers of t as like it 
would be in the limit Ay — > 0). Let us remind that the effective action is defined up to total derivatives, which can 
be omitted. Let us also stress that the above expression is not a proper-time expansion, since the coefficients Ciii 2 ...i n 
are functions of t. 

The integrals over the proper-time t can be simply evaluated and reduced to combinations of some set of elementary 
integrals J n (m 2 ) 



dt 



r p(*,A 2 ), 



(11) 



where n is integer. However in some sense the result of such integrations does not mean much by itself. One has 
to put the series in a form in which every term of the asymptotic expansion is chiral invariant if W[Y] does. This 
natural requirement selects from the infinite number of possibilities for rearrangements in Eq.([To|) only one as the 
correct answer. The main problem here is to find the algorithm which automatically yields a chiral invariant grouping 
for the background fields as well as the mass dependent factors before them. This problem has been solved recently 
on the basis of the novel algorithm for resummations inside the starting expansion (pp|). The resummations are 
determined by the recursion formulas 
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J t (m 2 ) - J,(m?) = ]T [J* + „(mf) - (-1)" J l+ ^(m?)] . (12) 
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In the case under consideration one has to factorize the mass dependent factors ij, manipulating formula (12). Let 



us stress that only the following choice for Ii leads to invariant coefficients at in the resulting asymptotic series 

W ^= / ^E^-*^)- /^-5>K 2 ). (13) 
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The first four coefficients are equal to 



oo = l, a 1 = -Y, a 2 = ^- + ^A 3 r+^=(Ai3 + A 23 )A 8 r, 

Y 3 1 1 
a 3 = — gj- - ^ A i2 ( A 3i + A 32 ) A 3 F + — y= [A 13 (A 2 i + A 23 ) + A 23 (Ai 2 + A i3 )] X 8 Y 

+ -L (A 31 + A 32 ) X 8 Y 2 + iA 21 A 3 F 2 - ^(dY) 2 . (14) 

To obtain this result we used relations between Ei matrices and U(3) hermitian generators, Ao,A 3 and As- One can 
consider the resulting series an inverse mass expansion, since Ii+i ~ m ^ 21 f° r £ > 1- However one should remember 
that the asymptotic coefficients a* depend on mass differences. This dependence is completely fixed by the symmetry 
requirements. Indeed, it can be verified that if the operator D^D is defined to transform in the adjoint representation 
5(D^D) = i[uj, D], the coefficient functions a,; are invariant under the global infinitesimal chiral transformations 
with parameters u) = a + 75/?. The present result is in agreement with the standard Schwingcr-DcWitt inverse mass 
expansion, when mi = m 2 = m 3 . For the case with Ay 7^ our formula ( p"3| ) is a new asymptotic series which can 
be used to construct the low-energy EFT action when the local vertices are induced by one-loop diagrams involving 
heavy particles with different masses. In the low energy QCD, for instance, our approach can be used for the Nambu 
- Jona-Lasinio type models with SU(3) x SU(3) chiral symmetry fll5|j , or for heat kernel one-loop renormalizations 
in the framework of the chiral expansion || . The difference in the masses of the nonstrange and strange constituent 
quarks is large enough, m u , ~ 330 MeV and m s ~ 510 MeV, so that one may expect essential numerical deviations 
from the results obtained by other methods. Finally the method can also be generalized to accommodate minimal 
couplings of the loop particles to gauge fields. Work along these lines is in progress. 
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